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基于 严格 极限 平衡 法 的 边 坡 安全 系数 上 下 限 解 
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摘 要 :严格 极限 平衡 法 满足 所 有 的 平衡 条 件 , 通 常 计 算 精 度 较 高 ,但 不 能 保证 所 有 条 块 界面 上 均 
不 违背 摩尔 -库仑 准则 ,也 不 能 保证 充分 发 挥 其 抗 剪 强度 ,因此 在 理论 上 既 不 是 下 限 解 .也 不 是 上 限 
解 , 且 计算 结果 一 定 程 度 上 依赖 于 条 块 间作 用 力 关系 的 假设 。 将 安全 系数 作为 设 定 值 , 以 滑 面 正 应 
F 力 为 未 知 变量 ,将 水 平地 震 力 影响 系数 作为 目标 函数 ,根据 边 坡 整 体 的 水 平 力 平衡 条 件 , 建 立 目 标 
函数 表达 式 ; 再 根据 边 坡 局 部 力 平 衡 条 件 及 条 块 间接 触 面 .滑动 面 上 的 破坏 准则 建立 不 等 式 约束 条 
件 ,根据 边 坡 整体 的 坚 向 力 与 力 珑 平衡 条 件 建立 等 式 约束 。 目 标 函 数 及 约束 条 件 构成 标准 线性 规 
划 数 学 模型 ,应 用 单纯 形 法 可 以 精确 得 到 设 定安 全 系数 下 的 使 边 坡 达到 极限 平衡 状态 的 最 大 与 最 
OO 小 水 平地 震 力 影响 系数 。 通 过 和 迭代 求解 ,得 到 既定 水 平地 震 力 影响 系数 下 的 基于 严格 极限 平衡 的 
边 坡 安 全 系数 上 下 限 解 答 。 算 例 计算 结果 表明 ,这 种 安全 系数 上 下 限 解 范 围 一 般 小 于 5% , 且 条 块 
间作 用 力 不 违 背 破坏 准则 且 充 分 发 挥 抗 剪 强度 ,因此 计算 结果 更 可 靠 , 且 理论 上 更 为 严密 。 
关键 词 :严格 极限 平衡 ;安全 系数 ;水 平地 震 力 影响 系数 ;上 下 限 解 ; 线 性 规划 ;优化 
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The upper and lower-bound solutions of slope safety factor 
based on rigorous limit equilibrium method 
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Abstract :The rigorous limit equilibrium method satisfies all of the equilibrium conditions ,commonly with 
high accuracy. However,such a solution may represent neither the lower solution ,nor the upper one ,as the 
forces acting upon the interslice surfaces with some assumptions cannot ensure non-violation of failure cri- 
terion and full mobilization of those shear strengths ,the results of which ,to some degree ,are dependent up- 
on the assumptions made about the relationships between the interslice forces. In the paper ,with the value 


of the factor of safety initially assumed ,the normal stresses acting on the slip surface are taken as the un- 
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knowns and the coefficient of horizontal seismic forces as the objective function ,the equation of the objec- 
tive function is established according to the horizontal force equilibrium for the whole body of sliding mass. 
Then ,the inequality constraint conditions are established according to the equilibrium condition of local 
sliding body and the failure criteria for both interslice surfaces and the slip surfaces ,and the equality con- 
straint condition is further established according to the vertical force and moment equilibrium conditions for 
the whole body of sliding mass. The equation of objective function as well as these constraint conditions 
thus constitutes a standard linear programming program ;the maximum and minimal values of the coefficient 
of horizontal seismic forces, causing the sliding mass into rigorous limit equilibrium condition ,can be pre- 
cisely determined by using the simplex method. The final upper- and lower bound solutions for values of 
the factor of safety associated the prescribed coefficient of horizontal seismic force can be obtained by a- 
dopting an iterative procedure. The results of example study show that the difference between the upper 
and lower bounds of the factor of safety is within the range of 5% of the average value. Since the interslice 
forces on the all the interslice surfaces do not violate the failure criterion and,the shear strength of which 
js ,to the largest extent, mobilized. Therefore , the results of present method are reasonable and theoretically 
~ rigorous. 


~ Key words: rigorous limit equilibrium; safety factor; upper-bound solution; lower-bound solution ; linear 


= programming ; optimization 


中 〇 尽管 有 限 元 等 数值 方法 在 岩 土工 程 领域 得 到 广 
泛 应 用 ” ,但 极限 平衡 法 仍 是 边 坡 工程 中 最 常用 分 
桥 尖 法 ,并 被 各 类 工程 勘察 设计 规范 纳入 必 选 方 
潜 S)  。 早 期 极限 平衡 法 如 瑞典 条 分 法 '9 .简化 Jan 
b 全 法 美国 陆军 工程 师 团 法 已 等 只 满足 部 分 力 或 
力 算 平衡 条 件 ,计算 结果 精度 较 差 (相对 于 满足 所 有 
平衡 条 件 的 严格 极限 平衡 法 ) 。 后 来 发 展 了 严格 极 
限 平 衡 法 如 严格 Janbu 法 0 、 Morgenstern-Price 
法 EE .Spencer 法 "及 基于 滑 面 正 应 力 修正 的 安 
全 又 数 显 式 解 法 ,这 些 严格 法 安全 系数 计算 精度 
大 的 提高 。 


力 假设 ,而 是 基于 滑 面 正 应 力 修正 的 极限 平衡 法 基 
本 思想 “ ,将 滑 面 正 应 力作 为 未 知 变 量 ,将 设 定 的 
安全 系数 下 使 边 坡 达 到 极限 平衡 状态 的 水 平地 震 力 
影响 系数 作为 目标 函数 ,根据 边 坡 整体 与 局 部 平衡 
条 件 及 滑 面 与 条 块 界 面 上 摩尔 库仑 破坏 准则 ,建立 
目标 函数 的 线性 表达 式 , 及 其 相关 线性 等 式 与 不 等 
式 约束 条 件 , 构 成 标准 的 线性 规划 问题 ,用 单纯 形 
法 ' ”求解 最 大 与 最 小 的 水 平地 震 力 影响 系数 。 通 
过 适当 迭代 方法 ,进而 求 出 基于 严格 极限 平衡 法 的 
边 坡 安全 系数 上 下 限 解 。 


但 必须 指出 ,尽管 严格 极限 平衡 法 满足 所 有 平 
衡 条 件 , 但 是 计算 结果 一 定 程度 上 还 依赖 于 条 间 力 
作用 方式 ,它们 计算 精度 高 是 基于 大 量 实际 计算 结 
果 经 验 总 结 得 出 的 ,没有 得 到 严格 理论 证 明 。 同 时 ， 
严格 极限 平衡 法 不 能 保证 所 有 条 间 力 都 能 不 超过 摩 
尔 库 仑 准则 所 限定 的 剪 切 力 极 限 值 ,因而 不 是 下 限 
解 ; 男 外 ,严格 极限 平衡 法 条 间 力 假设 不 能 充分 发 挥 
条 间 的 抗 剪 强度 ,因此 它 也 不 是 上 限 解 。 对 于 非 贺 
滑 面 不 同 严格 法 安全 系数 相差 也 会 达到 11% ” , 真 
正 的 严格 极限 平衡 法 安全 系数 合理 范围 至 今 仍 未 获 
得 ,因此 基于 严格 极限 法 的 上 下 限 解 的 研究 具有 重 
要 理论 意义 ,同时 缩小 严格 极限 法 安全 系数 可 能 范 
围 又 具有 很 大 的 工程 应 用 价值 。 

本 研究 不 采用 传统 的 严格 极限 平衡 法 中 的 条 间 


1 严格 极限 平衡 法 基本 方程 


对 于 给 定 滑 面 位 置 的 边 坡 (如 图 1 所 示 ) ,作用 
于 滑 体 上 的 力 有 滑 体 总 自重 WW 和 总 水 平地 震 力 
K.W , K. 为 水 平地 震 力 影响 系数 。 将 滑 体 划 分 为 n 
个 条 块 , 条 块 底面 宽度 为 b, (条 块 自 右 向 左 编号 )， 
则 条 块 底面 上 总 正 应 力 @ = [ol oa … 


ay,] 前 


? 


应 力 r = [ri mm … Ti ;及 和 孔隙 水 压力 w = 
[wu ws… wu,] 。 滑 面 上 各 条 块 底 面 的 有 效 内 摩 


氛 角 和 凝聚 力 分 别 为 p; 和 , 边 坡 安 全 系数 设 定 为 
有 ,由 摩尔 -库仑 强度 破坏 准则 得 
(0,—u,)tang,;, +c， 


二 辣 (1) 
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b.(1 + tanpitanas =W-— 1 S 2iCitana 
| i F. i i i F. A 二 i 

(3b) 


n 


| (1 十 Ftangitanas — x.)+ 


间 


i=1 


(tana, 一 Feang ,cr -7)|c， = 


F 2 Cy 7) tonal sw, ~ %.)] (3c) 


取 图 1 中 第 1 ~j 个 条 块 整体 (局 部 滑 体 ) 作为 
研究 对 象 ,如 图 2 所 示 , WW 为 第 1 ~j 个 条 块 整体 自 
村 重 , P; 和 7 分 别 为 作用 于 第 j 和 j+1 个 条 块 之 间接 
< 人。 触 面 上 的 法 向 力 和 切 向 力 , U; 为 作用 于 接触 面 上 的 
图 1 边 坡 滑 体 及 受 力 示意 图 水 压力 ,j = 1,2,… 和 -1。 


Fig.1 Diagram of sliding body and forces 

-严格 极限 平衡 法 要 求 滑 体 满足 所 有 力 和 力矩 平 
衡 乏 件 ， 对 于 二 维 边 坡 ,分 别 是 水 平 力 、 垂 直 力 和 力 
答 王 衡 。 根 据 滑 体 的 受 力 情况 ,对 滑 体 大 致 中 心 , 坐 
枯 Ox.,)) , 取 和 矩 , 有 


并 reore 一 并 visina = 天 于 


了 risina + wiieosa =W 图 2 局 部 滑 体 条 抉 上 的 作用 力 示意 图 
全 Fig.2 Forces on a local sliding body 
YIZ | 3 
ER 根据 水 平和 坚 直 方向 力 的 平衡 条 件 可 得 
n P., 二 KW + > b(tana, 一 Ftangi)o0 人 
py (Coilisinai - Tilicosai)Cyi — y.) = 0 ; Ee ” 
. De (4a) 
其 惠 : a ; 为 条 块 i 的 底面 倾角 ; Cxi ,yi) 为 第 i 个 条 Ff 和 
J 
二 而 中 心 全 标 ;条 据 底 而 长 度 为 = b/eose ,将 pi 
1 = bi/eosas 代入 以 上 各 式 , 则 整理 可 得 3 个 平衡 方 Ea | 
上 1 < 
程 , 印 赤 2 b,C,tana, (4b) 


3 > bi(7, 一 Cjitana',) 一 0 (2a) 
i=1 


2 线性 规划 问题 的 构建 


W-— > tet +0,)=0 (2b) 
n i 可 以 写成 


> (7,tana, + 0;) (Xx; — x )+ 


一 K, ol(} Fp -tanajui + 时 Yoo 
2 Caitanas -TD)Cxr -y)=0 (2c) (5) 
令 Ci= cwtang; ,将 式 (1) 代 入 式 (2) 的 3 个 由 式 (5) 可 以 看 出 , 边 坡 的 水 平地 震 力 影响 系 
平衡 方程 中 , 则 有 数 K, 与 滑 面 上 的 正 应 力 er 时 线性 关系 。 
>6(} np -aaa je = = KW -二 > C， 假设 边 坡 滑 面 固定 ,初始 设 定安 全 系数 F, 值 ， 


采用 线性 规划 方法 ,以 水 平地 震 力 影响 系数 K, 作为 
(3a) 目标 函数 ,以 滑 面 正 应 力 o 为 未 知 变量 。 以 天 . 为 目 
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衡 法 的 边 坡 安全 系数 上 下 限 解 InaXIV 王 1 


标 函 数 时 ,求解 得 到 的 天. 愈 大 ,说 明 奉 边 坡 达到 极 
限 平 衡 状态 时 所 需 的 水 平地 震 力 愈 大 ,对 应 的 安全 
系数 ,也 愈 大 ;反之 得 到 的 K, 愈 小 ,对 应 的 安全 系 
数 FF 也 愈 小 。 由 此 可 将 求解 边 坡 安全 系数 的 上 、 下 
限 解 问题 转化 成 求解 边 坡 处 于 极限 平衡 状态 时 水 平 
地 震 力 影响 系数 的 最 大 值 和 最 小 值 问 题 。 同 时 ,将 
滑 体 整体 满足 力 和 力矩 平衡 条 件 , 条 块 间接 触 面 不 
违背 摩尔 -库仑 准则 ,不 产生 拉 应 力 , 滑 体 满 足 运动 
许可 等 作为 约束 条 件 , 构 建 线性 规划 数学 模型 ,求解 
得 到 水 平地 震 力 影响 系数 的 最 大 值 和 最 小 值 ,进而 
得 到 基于 严格 极限 平衡 法 的 边 坡 安全 系数 上 限 解 和 
下 限 解 。 


2.1 目标 函数 建立 

王 式 (5) 中 , 令 

< A = [5 (Ftane, -tana | 

OO 1 

Sy b, (Flangp 一 tana | ee 
b, 人 tang, tana, ]] 要 

人 

CN ”FWE 

中 则 目标 函数 可 写成 

和 max K.=A'o+kK,, 
min 天 =Ar+K (6) 


2 723 线性 规划 问题 应 满足 的 约束 条 件 


2: 舌 1 各 条 间接 触 面 上 法 向 力 与 切 向 力 的 关系 不 
应 污 背 摩尔 -库仑 破坏 准则 
” 滑 体 作 为 整体 滑 移 ,不 应 在 条 块 间接 触 面 上 先 
发 生 破 坏 ,因此 作用 于 第 7 个 和 第 j+1 个 条 块 间接 
触 面 上 的 法 向 力 P; 和 切 向 力也 的 力学 关系 应 当 不 
违背 摩尔 -库仑 准则 , 即 有 
(Be U)tang; teiL, 
也 
其 中 :po 和 分 别 为 条 块 线 上 土 体 的 平均 内 摩擦 角 
和 黏 聚 力 ;L 为 作用 于 接触 面 上 的 水 压力 ;万 为 该 
条 块 线 的 长 度 。 
令 式 (7) 不 等 式 右边 为 卫 , 即 
(Pi;- UV)tangp， +toL, 
2 
则 DP 表示 在 临界 状态 下 ,由 摩尔 -库仑 准则 求 
得 的 条 间接 触 面 上 切 向 力 的 理论 值 ,根据 上 述 约束 
条 件 ,条 间接 触 面 上 切 向 力 的 实际 值 7 应 不 大 于 理 


Ts ! 


(J=1,2,…,n-1) (7) 


T= 
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论 值 五 , 即 式 (7) 可 表示 为 7<7。 
由 式 (4a) (4b) 和 式 (6) (7) 整 理 可 得 
(WA -6 ) os<K -Kon, (8) 
式 中 ， 
-Wtangi 
— W,tang, 
i= 2 ， 
Wang, 
Fen 0 0 “0 
2 Cy 0 “… 0 


SS 
1 


Ci 3 a … 0 


1 1 
一 全 人 68 2 人 一 


1 
Cy = bl F.+tang,tana, +tanpitanay 一 Fanpjtang 


(=1.2 ,jj=12 ,1-1), 
a bt a 一 人 
FW + -0, + b,C,tana, 
， tang’; 2 2 
Kk = 一 下 用 +c',L, Tp > DC 十 > b,Ctana, 
S sl 站 二 下 
t 2 
-FW+e'L - 人 YbC + YbiCitana, 
S a Ss 一 


若 需 所 有 条 块 间 接触 面 上 法 向 力 和 切 向 力 的 关 
系 均 满 足 该 约束 条 件 , 则 式 (8) 共计 由 守 - 1 个 方程 
组 成 。 
2.2.2 条 块 间 不 产生 拉 应 力 

土 体 抗 拉 强度 很 小 ,因此 条 块 间 不 应 产生 拉 应 
力 , 对 于 滑 面 人 入 口 可 能 出 现 的 土 体 拉 应 力 情况 ,可 设 
置 张 拉 裂缝 予以 消除 。 以 图 2 所 示 忆 的 方向 为 正 
方向 , 则 应 有 条 间接 触 面 上 的 有 效法 向 推力 已 =- 
U > 0 ,由 式 (4a) 和 式 (6) 整 理 可 得 

- (WA +AVo <K, (9) 


其 中 ， 


于 =[W， 有 一 ] 


J 


| (: t | 
1| tanal 一 到 tanpl 
Feang) 


1 
| 2 ( tana， 一 产 ee . 


1 
bl tana, 一 amp] 


J 
K, 二 天 全 ) 过 二 DC a U; » 


CI=lo 0o,; oj (j=1,2,.…,n-1) 
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式 (9) 同 样 由 nn-1 个 方程 组 成 。 来 求解 该 线性 规划 问题 。 

2.2.3 滑 体 满足 运动 许可 求解 水 平地 震 力 影响 系数 K. 最 大 值 ,首先 假设 
滑 体 及 条 块 间 接触 面 上 的 作用 力 需 满足 运动 许 。” 初始 安全 系数 ,由 上 述 构建 的 目标 函数 方程 (6)， 


可 , 即 如 图 2 所 示 的 7 方向 为 正方 向 ,还 应 有 7 > 
0, 由 式 (4b) 整 理 可 得 


Aro < 天 7 (10) 
其 中 ， 
oli + Flanputana | 
b (! de tana ) 
Ns 2 a Pp) 2 ， 
1 
| to + Flangpilanes 
, Te 
天 7 = W, = 天 之 b;C.tana, ， 
sr=[c o, oj] (j=1,2,…,n-1) 
式 如 0) 由 n -1 个 方程 组 成 。 
滑 体 满足 坚 直 力 平衡 条 件 
根据 滑 体 竖 直 力 的 平衡 方程 , 邻 
r 1 
bl(1 十 Fanptanan ) 
有 Po b,(1 十 Ftangatanae) 
b (1 + ne tana ) 
L 天 下 n n 
E K,=W-— 元 > biCitana, 
& Fr i=1 
己 则 式 (3b) 可 记 为 
KL, lg = = kK, (11) 
2.2. 5 滑 体 满 足 力 给 平衡 条 件 


恨 据 滑 休整 休 的 力 和 平和 方程, 仿 令 
J = 二 [ui, 2, ; HA] 


Li tana; | (x; —x.) + 
Ff op; oj] Xi — Ne 
1 
(bane -Ftangi) (9 =7) (i=1,2,.,n), 
K i 一 )— tana,.(x —x)] 
M Fr i=1 : yi Ye 了 i 他 


则 式 (3e) 可 记 为 


Lno = Ky (12) 


3 求解 步骤 


根据 上 述 构建 的 数学 模型 ,可 编写 计算 机 程序 


投稿 网 站 :http://cjam. xjtu. edu. cn 


以 及 滑 体 和 条 块 间接 触 面 应 满足 的 各 项 约束 条 件 方 
程式 (8) ~ (12)。 应 用 单纯 形 法 ,可 以 求解 得 
到 天. 。 通 过 对 比 求解 得 到 的 天. 与 边 坡 原本 给 定 的 
天 6 ,调整 程序 中 代入 的 初始 安全 系数 , 若 K。 < 
Ko , 则 说 明 假 定 的 初始 安全 系数 FR, 比 实际 的 安全 
系数 大 ,为 使 大. 增 大 , 则 需 减 小 代入 的 安全 系数 , 反 
之 同 理 。 知 干 次 迭代 后 , 当 |K. - K,| 二 8 时 ,说 明 
KK, 与 Ko 充分 接近 , 则 此 时 的 了 . 即 为 解答 。e 视 精度 
要 求 决 定 , 对 于 一 般 工 程 问 题 ,可 取 为 0. 005 ~0.01。 
由 于 水 平地 震 力 影响 系数 K, 与 滑 面 上 的 正 应 力 er 
呈 线 性 关系 ,上 述 求解 过 程 往往 只 需要 迭代 几 次 便 
可 达到 精度 要 求 。 求 解 水 平地 震 力 影响 系数 天 最 

小 值 过 程 与 之 相同 ,此 处 不 再 歼 述 。 


4 算 例 分析 


以 下 通过 两 个 算 例 的 求解 ,并 与 传统 极限 平衡 法 
的 计算 结果 作对 比分 析 ,验证 本 研究 方法 的 合理 性 。 


4.1 算 例 1 


如 图 3 所 示 的 边 坡 算 例 辣 , 边 坡 剖面 、 滑 动 面 
及 土 体 各 项 参数 标 于 图 中 。 本 研究 计算 分 析 如 图 中 
所 示 的 两 个 滑 面 , 滑 面 1 为 圆 弧 滑 面 , 滑 面 2 为 非 圆 
弧 滑 面 ,两 个 滑 面 均 为 固定 滑 面 ,是 与 最 小 安全 系数 
无 关 。 在 本 算 例 中 ,既定 的 水 平地 震 力 系数 KK。= 
0 ,分 别 采用 传统 算法 、 其 他 学 者 经 典 算法 '“ 和 本 方 
法 计算 所 得 的 安全 系数 见 表 1。 以 滑 面 1 的 计算 结 
果 为 例 ,用 本 方法 计算 所 得 的 滑 面 正 应 力 大 小 、 各 条 
块 间接 触 面 上 的 推力 大 小 以 及 屈服 条 件 与 简化 的 
Bishop 法 的 计算 结果 对 比分 析 分 别 见 图 4 ~6。 


R=24.384 m 


y=18.852 kN/m’ 
C=28.74 kPa 


涓 面 2， 非 圆 弧 


图 3 算 例 1 示意 图 


Fig.3 No.1 example slope and parameters 
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表 1 算 例 的 计算 结果 与 比较 
Tab.1 Calculated results and comparison 


本 


序号 计算 方法 
滑 面 1, 圆 弧 滑 面 2, 非 圆 弧 

1 简化 Bishop 法 2.080 国 
2 Spencer 法 2.073 1.373 
3 Morgenstern & Price 法 2.076 1.370 
4 ”基于 滑 面 正 应 力 修正 法 。 2.061 1.392 
5 ”本 研究 方法 (下 限 解 ) 2.062 1.371 
6 。 本 研究 方法 (上 限 解 ) 2.082 1.395 


~ 


上 “本 研究 方法 (下 限 解 bp 
$ 本 研究 方法 (上 限 解 ) 、、 


0 : 
-5.760 4.584 28.720 


14.928 
xX/m 


21.824 


图 4 条 块 间接 触 面 推力 大 小 对 比 


Fig.4 Comparison of the thrust lines 


简化 Bishop 法 盖 了 


0 4 一 
-760 2.860 13.204 23.548 28.720 


xX/m 


图 5 滑 面 上 正 应 力 分 布 
ig.5 Normal stress distribution on sliding surface 
250 


200 


日 
F150 
& 100 1 
BB ” 实际 值 7 

50 

j 

-5.760 4.584 14.928 21.824 28.720 

xX/m 
(a) AF. 下 限 解 时 

250 

200 理论 值 7' 
150 
又 100 5 
S 实际 值 7、 


50 


21.824 


0 
-5.760 4.584 14.928 


xX/m 
(b) 上 限 解 时 
图 6 条 间 切 向 力 的 理论 值 与 实际 值 对 比 


Fig.6 Comparison between the theoretical value and 


28.720 


actual value of the tangential force on interslice surfaces 
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由 表 1 可 见 , 对 于 圆 弧 滑 面 ,本 人 研究 方法 计算 的 
安全 系数 的 下 限 解 与 基于 滑 面 正 应 力 修正 法 所 得 结 
果 十 分 接近 ,其 他 严格 条 分 法 计算 结果 均 处 于 本 研 
究 方 法 上 下 限 解 的 范围 内 ;对 于 非 贺 孤 滑 面 ,本 研究 
方法 计算 的 安全 系数 的 下 限 解 与 Morgenstern & 
Price 法 (以 下 简称 “M-P 法 ”) 所 得 结果 十 分 接近 ,其 
他 严格 条 分 法 计算 结果 也 均 处 于 本 研究 方法 上 下 限 
解 的 范围 内 。 从 图 6 中 可 以 看 出 , f, 取 上 下 限 解答 
时 ,条 块 间接 触 面 上 的 受 力 情况 均 不 违背 摩尔 -库仑 
准则 。 由 本 研究 方法 通过 线性 规划 得 到 的 安全 系数 
上 限 解 和 下 限 解 相差 约 1% 。 


4.2 算 例 2 


如 图 7 所 示 的 边 坡 算 例 ,图 中 一 般 滑 面 和 圆 弧 
滑 面 均 为 固定 请 面 , 且 与 最 小 安全 系数 无 关 , 各 土 层 
的 参数 属性 见 表 2。 分别 采 用 简化 的 Bishop 法 , 改 
进 的 M-P 法 "和 本 研究 方法 求解 。3 种 方法 计算 
所 得 结果 对 比 见 表 3。M-P 法 的 基本 假设 是 条 块 间 
法 向 力 与 前 切 力 的 比值 用 条 间 力 函数 f(x) 和 待定 
系数 和 的 乘积 表示 。 在 改进 的 M-P 法 中 选择 的 条 
间 力 函数 


f(x) = sin’ | mr (7 = 2) ] (14) 


b-a 


式 中 :a 和 4 分 别 为 边 坡 的 左右 边界 ;yw 和 vw 为 非 负 
参数 ,一 般 可 取 =0 ~ 5.0 ,> = 0.5 ~ 2.0。 采用 
不 同 参数 求解 得 到 的 待定 系数 A 和 边 坡 安全 系数 
下 如 表 4 所 示 ( 以 一 般 滑 面 , 湿 坡 ,K。 =0. 1 为 例 ) 。 


30.0 


图 7 算 例 2 示意 攻 
Fig.7 No.2 example slope 
表 2 边 坡 土 层 参数 属性 


Tab.2 Parameters of every stratum 


土 层 y/(kKN.m’) C/ kPa p/(°) 
OD 18.2 20.0 32.0 
@) 18.0 25.0 30.0 
®@ 18.5 40.0 18.0 
@ 18.8 40.0 28.0 
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表 3 几 种 方法 计算 的 安全 系数 比较 


Tab.3 Comparison of calculated results 


圆 弧 滑 面 一 坡 滑 面 


干 坡 湿 坡 干 坡 湿 坡 
Ko 0 0.1 0 0.1 0 0.1 0 0.1 
简化 的 Bishop 法 2.618 1.857 1.681 1.181 一 一 二 
M-P 法 * 2.618 1.969 1.752 1.315 2.197 1.598 1.615 L171 
本 研究 方法 (下 限 解 ) 2.553 1.937 1.707 1.270 2. 156 1.415 1.462 1.038 
本 研究 方法 (上 限 解 ) 2.619 1.981 1.755 1.330 2.313 1.654 1.701 1.235 
注 :“ 本 表格 中 M-P 法 的 计算 结果 为 条 间 力 函数 是 常数 时 。 
通过 表 3 中 的 计算 结果 对 比 发 现 ,大 部 分 简化 一 Morgenstern & Priced 


. 本 研究 方法 (下 限 解 ) 
一 本 研究 方法 (上 限 解 ) 2 区 


的 Bishop 法 和 条 间 力 函数 为 常数 时 的 M-P 法 计算 
所 得 的 安全 系数 均 处 于 本 研究 方法 上 下 限 解 的 区 间 
内 ; 圆 弧 滑 面 干 坡 情况 下 KK =0.1 时 ,简化 Bishop 
法 的 计算 结果 小 于 本 研究 方法 下 限 解答 ,两 者 相差 J 
4 本 5 ; 圆 弧 滑 面 湿 坡 情况 下 简化 Bishop 法 计算 结 A 


P/(KN-m') 


者 用 于 本 研究 下 限 解 答 ， 当天 o 按 表 中 取 值 时 ,分 别 图 9 一 般 滑 面 三 种 方法 推力 大 小 对 比 ( 湿 坡 ,Ko =0.1) 
= 1.2% 和 7.3% ;本 人 研究 方法 的 下 限 解 通常 都 小 Fig.9 Comparison of thrust lines on interslice surfaces 
其 余 余 方 法 的 计算 结果 , 即 利用 本 研究 方法 可 以 找 for three methods (wet slope,K.o =0.1) 
人 表 4 M-P 法 不 同 条 间 力 函数 选取 时 的 安全 系数 
以 一 般 滑 面 湿 坡 条 件 下 Ks =0.1 这 一 工 况 为 (以 一 般 滑 面 , 湿 坡 ,K.。=0.1 为 例 ) 


Tab.4 Safety factors by M-P method 


人 对 比 表 3 和 表 4 中 的 计算 结果 可 以 发 现 , 当 M-P 
法 阅 条 和 间 力 函数 取 不 同 参数 时 ,计算 结果 均 处 于 本 


( general slip surface,wet slope, Ku =0.1) 


研究 方法 上 下 限 解 的 区 间 内 。 4 2 5 
一 3 种 方法 计算 所 得 的 正 应 力 分 布 情况 见 图 8, 推 .3 ed el 

湖人 小 对 此 兄 图 9, 采 用 本 研究 方法 条 块 线 上 法 向 了 3 0 en 

力 独 切 向 力 满足 摩尔 - 库伦 准则 情况 见 图 10。 L :909 Ho 
. 一 Morgenstern & Price 法 2 0.3671 1.152 
本 研究 方法 (下 限 解 ) 入 1 0.5 0.3615 1.165 


一 -本 研究 方法 (上 限 解 ) 


1 0.3550 1.137 

2 0.446 1 1.156 

2 0.5 0.4461 1.162 

1 0.4332 1.110 

图 8 一 般 滑 面 3 种 方法 正 应 力 分 布 对 比 ( 湿 坡 ,K.。 =0.1) ee 
Fig.8 Comparison of normal stress distribution on slip surface ” 3 人 0 
for three methods (wet slope,K,, =0.1) 0 Ee 

此 外 ,值得 注意 的 是 ,在 不 同 工 况 下 ,本 人 研究 方 Vb ee 

法 计算 所 得 的 上 下 限 解 答 之 间 的 差异 也 不 尽 相 同 ， 4 0 Od Eo 
二 者 相差 最 小 仅 为 2% ,最 大 约 17% , 且 如 图 8 中 所 , A 0 
示 , 正 应 力 分 布 波动 较 大 ,存在 多 个 突变 点 ,这 是 由 g 0 been 
于 本 研究 中 的 滑 面 正 应 力 er 是 离散 的 ,并 未 对 其 进 > 3 CI wo 
行 其 他 条 件 的 约束 , 故 本 研究 方法 是 满足 数学 解 和 下 
答 的 。 2 0.8209 1.219 
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EF 衡 法 的 边 坡 安全 系数 上 下 眼角 InaxX | Ve 作 期 |- 


理论 值 7 


0 


7kN m 


X/m 


(a) 下 限 解 时 


理论 值 7' | 


TKN* m') 


实际 值 7 了 } 


X/m 


(b) 已 上 限 解 时 


一 般 滑 面条 间 切 向 力 理论 值 与 实际 值 对 比 
( 湿 坡 ,K.,=0.1) 


g.10 Comparison between the theoretical and actual 


ue of the tangential force on interslice surfaces of slope 
with a general slip surface (wet slope ,KK =0.1) 

C9 车 结合 实际 情况 ,可 同时 将 滑 面 正 应 力 分 布 曲 
线 为 凸 曲线 作为 一 组 约束 条 件 , 则 所 得 上 下 限 解答 
的 珊 间 范围 会 大 大 缩小 。 仍 以 算 例 2 中 的 一 般 滑 面 
湿 披 条 件 下 K。 = 0. 1 为 例 ,未 对 正 应 力 分 布 约束 
-起 下 限 解 答 分 别 为 1. 235 和 1. 038, 相差 约 
% ;对 正 应 力 分 布 加 以 约束 时 ,上 下 限 解答 分 别 
097 和 1. 056 ,相差 约 4% 。 

-人 此 时 的 正 应 力 分 布 曲线 及 摩尔 -库伦 准则 遵循 
情况 与 M-P 法 的 对 比 可 分 别 见 图 11 和 图 12 ,可 以 
看 出 ,M-P 法 的 正 应 力 分 布 存在 “凹陷 ”的 情况 , 且 当 
nw 和 vw 取 值 较 大 时 ,条 间接 触 面 上 的 法 向 力 与 切 向 
力 并 不 能 完全 满足 摩尔 -库仑 准则 ,因此 可 以 认为 此 
时 的 安全 系数 也 是 不 合理 的 。 


Morgenstern & Price 法 
ee 本 研究 方法 (下 限 解 ) 
-一 -本 研究 方法 (上 限 解 ) 二 <、、 


o/kPa 


xX/m 


约束 后 的 滑 面 正 应 力 分 布 情况 
(一 般 滑 面 , 湿 坡 ,K.。=0.1) 


图 11 


Fig.11 Normal stress distribution on slip surface after 


constraint (general slip surface ,wet slope,K.o =0.1) 


刊 ss7 


实际 值 7 


“、、、 ”理论 值 7 
>< 


TKN* nm 


X/m 
M-P 1=5, v=2 


图 12 M-P 法 条 间 力 函数 j=5,v=2 时 
条 间 切 向 力 理论 值 与 实际 值 对 比 


Fig. 12 Comparison between the theoretical and actual 


value of the tangential force on interslice surfaces 


by M-P method (1 =5,v=2) 


尽管 大 量 计算 结果 表明 不 同 严 格 极限 平衡 法 计 
算 结 果 在 工程 意义 上 是 非常 接近 的 ,通常 不 超过 
10% ,但 至 今 仍 无 法 得 到 严格 极限 平衡 法 的 上 下 限 
范围 。 无 论 是 探究 严格 极限 平衡 法 的 上 下 限 范围 ， 
还 是 在 保证 不 违背 破坏 准则 与 运动 许可 条 件 情况 
下 ,将 基于 严格 极限 平衡 法 的 边 坡 安全 系数 上 下 限 
范围 进一步 减 小 ,其 在 理论 与 工程 应 用 方面 都 具有 
重要 意义 。 本 研究 在 严格 极限 平衡 法 理论 框架 下 ， 
以 滑 面 正 应 力 为 未 知 变量 , 以 水 平地 震 力 为 目标 函 
数 , 根 据 滑 体力 与 力矩 平衡 条 件 以 及 滑 面 、 条 块 界面 
的 剪 切 强度 准则 不 抗 拉 条 件 和 运动 许可 机 制 ,建立 
标准 线性 规划 问题 ,利用 标准 的 单纯 形 算法 可 直接 
求解 边 坡 在 设 定 的 安全 系数 下 ,水 平地 震 力 影响 系 
数 最 大 与 最 小 值 。 应 用 简单 的 迭代 改进 方式 ,进而 
求 出 给 定 的 水 平地 震 力 影响 系数 下 的 边 坡 安全 系数 
上 下 限 解 。 通 过 算 例 的 计算 对 比分 析 , 得 到 了 以 下 
几 点 结论 。 

1) 将 求解 边 坡 安 全 系数 的 上 、 下 限 解 问题 转化 
为 求解 极限 平衡 状态 下 水 平地 震 力 影响 系数 的 最 大 
与 最 小 值 问题 的 好 处 是 ,将 复杂 非 线性 规划 问题 转 
化 为 唯一 确定 解 的 标准 线性 规划 问题 。 

2) 计 算 结 果 表明 ,对 给 定 的 圆 弧 与 非 贺 激 滑动 
面 都 可 以 得 到 满足 所 有 约束 条 件 的 安全 系数 上 下 限 
解 , 对 圆 孤 滑 面 ,上 下 限 解 相差 小 于 5% ; 对 于 非 圆 
弧 滑 面 , 上 下 限 解 最 大 相差 可 接近 20% 。 

3) 如 果 进 一 步 将 滑 面 正 应 力 分 布 形 状 进行 约 
束 , 使 其 不 会 出 现 突变 , 则 非 圆 弧 滑 面 的 安全 系数 上 
下 限 差异 也 能 减 小 到 5% 以内。 

4) 传 统 的 简化 Bishop 法 与 Morgenstern-Price 法 
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engineering manual [ R |]. Vicksburg, Mississippi: Army Corps of 


大 多 数 情 况 下 安全 系数 在 本 人 研究 方法 的 上 、 下 限 解 


之 内 ,但 有 时 也 在 此 范围 之 外 , 且 条 间 力 出 现 违 背 破 
坏 准 则 的 情况 ,也 不 能 从 理论 上 充分 发 挥 条 块 界 间 
的 抗 剪 强 度 。 

5 ) 本 研究 基于 严格 极限 平衡 法 , 既 保 证 了 竖 直 
条 块 间 的 作用 力 不 违 背 破 坏 准则 ,又 在 严格 极限 平 
衡 条 件 下 , 尽 可 能 发 挥 条 间 的 抗 剪 强度 ,因此 理论 上 
更 加 严密 ,计算 结果 更 为 可 靠 。 

本 研究 是 基于 竖 直 条 块 的 严格 极限 平衡 法 ,但 
不 能 保证 其 他 斜 条 块 界面 的 作用 力 严 格 遵守 破坏 准 
则 ,如 果 进 一 步 增加 所 有 可 能 斜 条 块 界面 的 约束 条 
件 ,有 可 能 将 现在 的 5% 上 下 限 解 范围 进一步 减少 。 
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